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ABSTRACT 

One of the conundrums in extragalactic astronomy is the discrepancy in ob- 
served metalhcity distribution functions (MDFs) between the two prime stellar 
components of early-type galaxies — globular clusters (GCs) and halo field stars. 
This is generally taken as evidence of highly decoupled evolutionary histories 
between GC systems and their parent galaxies. Here we show, however, that 
new developments in linking the observed GC colors to their intrinsic metallici- 
ties suggest nonlinear color-to-metallicity conversions, which translate observed 
color distributions into strongly-peaked, unimodal MDFs with broad metal-poor 
tails. Remarkably, the inferred GC MDFs are similar to the MDFs of resolved 
field stars in nearby elliptical galaxies and those produced by chemical evolution 
models of galaxies. The GC MDF shape, characterized by a sharp peak with 
a metal-poor tail, indicates a virtually continuous chemical enrichment with a 
relatively short timescale. The characteristic shape emerges across three orders 
of magnitude in the host galaxy mass, suggesting a universal process of chemical 
enrichment among various GC systems. Given that GCs are bluer than field 
stars within the same galaxy, it is plausible that the chemical enrichment pro- 
cesses of GCs ceased somewhat earlier than that of field stellar population, and 
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if so, GCs preferentially trace the major, vigorous mode of star formation events 
in galactic formation. We further suggest a possible systematic age difference 
among GC systems, in that the GC systems in more luminous galaxies are older. 
This is consistent with the downsizing paradigm whereby stars of brighter galax- 
ies, on average, formed earlier than those of dimmer galaxies; this additionally 
supports the similar nature shared by GCs and field stars. Although the sample 
used in this study (the HST ACS/WFC, WFPC2, and WFC3 photometry for 
the GC systems in the Virgo galaxy cluster) confines our discussion to R < Re 
for giant ellipticals and < 10 Re for normal ellipticals, our findings suggest that 
GC systems and their parent galaxies have shared a more common origin than 
previously thought, and hence greatly simplify theories of galaxy formation. 

Subject headings: galaxies: evolution — galaxies: elliptical and lenticular, cD — 
galaxies: individual (M84, M87, NGC 147, NGC 3377, NGC 3379, NGC 5128) 
— galaxies: star clusters — globular clusters: general 



INTRODUCTION 



Globular cluster (GC) systems are always present in large galaxies. It is generally 
believed that GCs form when starbursts occur in galaxies, and as 'fossil re cords', contain 



vital information on the formation and evolution of their parent galaxies (ISear 



1978 



Harrislll99l[ kshman fc zjll998[ IWest et al.l[20o3 ISrodie fc Stradeilbood : 



e &: Zinn 



Lee et al. 



2010al ). Despite their close interplay between GCs and field stars, the comparative studies 
have uncovered a fundamental difference in the observed shapes of their metallicity distribu- 
tion fui]£tion^_(MDF^^ galax- 
ies faarris fc Harrislbood : [lleikuba et aPboOsi [Harris et aLll2007al lbl: Isird eraPboioh . The 
cause of the discrepancy between these two pri me stellar conaponents of galaxies has been the 
topic of much interest both on theoretica l (e.g., Beaslev et al.ll2002 : Pipino. Puzia. fc Matteucci 



2007 ) and observatio nal grounds (e.g., iForbes fc Forte 



2001 



Forte. Faifer. fc Geislei 



2005 



20071 : iLiu et al.ll201l[ ) because the disagreement signifies highly decoupled evolutionary paths 
between GC systems and their parent galaxies. 

The most common technique for measuring metallicities of a substantially large sample 
of GCs is photometry - obtaining their broadband colors, although it is no substitute for 
spectroscopy. Because GCs in the Milky Way and other galaxies are usually older than 
10 Gyr and age does not strongly affect GC broadband colors of GCs this old, the main 
parameter governing GC colors is metallicity. Indeed, the overall, first-order feature of the 
color-metallicity relations (hereafter "CMRs") is that GC colors scale linearly with their 
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metallicities. Empirical relationship s between the most often used color , V — /, and [Fe/H], 



fitted mainly to the Galactic GCs ( 


Couture. Harris. & AUwrieht 


1990; 


Kissler-Patie et al. 


1997. 


1998b: 


Kundu & Whitmore 


1998; 


Barmbv et al. 


2000). are aDDroximatelv linear. Us- 


inp; the more metallicity-sensitive colors C — or C — R, both 


Harris & Harris 


(I2OO2 


) and 



Cohen. Blakeslee. fc Cotel ( l2003l ) found a mildly quadratic or broken linear relationship be 



tween [Fe/H] and color to be a better fit, thus improving on the relation of iGeisler fc Forte 



fll990[ ). 



observation of bimoda 


ity in GC color distributions (e.er 


.. Zepf & Ashman 


1993; Ostrov. Geisler. & Forte 


1993; 


Whitmore et al. 


1995 


Lee. Kim. & Geisler 


1998: 


Gebhardt & Kissler-Patie 1999: 


Harri 




2001; 


Kundu & Whitmore 


2001; 


Larsen et al. 


2C 


01; 


Pene et al. 2006; 


Harris et al. 


2006; 




Jordan et al. 2009; Sinnott et al. 


2010; Liu et al. 


201 


1) has been translated into bimodal- 





ity of their MDFs. This is where a sharp distinction between GCs and field stars takes 
place; independent studies via direct photometry of spatially-resolved constituent field stars 
in a dozen nearby galaxies have shown that their MDFs have, in general, strongly-peaked, 
unimodal [Fe/H] distributions with broad metal-poor tails. 

More recent observations and modeling of old star clusters, however, suggest that the 
relations between metallicity and broadband color for GCs have a subtle, second-order fea- 
ture; the y appear to be no nlinear with a quasi-inflection at intermediate metallicities. For 
instance, iPeng et al.l (|2006[ ) presented an empirical relationship between the g — z colors and 
spectroscopic metallicities for GCs in the Milky Way and the giant elliptical galaxies, M49 
and M87 (see their Figure 11). With this dataset, they showed that the relationship between 
[Fe/H] and g — z is steep for [Fe/H] < —0.8, shallow up to [Fe/H] ~ —0.5, and then possibly 
steep again at higher metallicities. 



Independently, lYoon. Yi. fc Led (120061 . hereafter Paper I) presented a theoretical metallicity- 
to-color relationship that has a significant inflection and thus reproduces well the observed 
feature. This nonlinear nature of the relation between intrinsic metallicity and its proxy, col- 
ors, may hold the key to understanding the color bimodality phenomenon. Paper I showed 
that the wavy feature projects equidistant metallicity intervals near the quasi-inflection point 
onto larger color intervals, and thus produces bimodal GC color distributions when the un- 
derlying distribution in [Fe/H] is broad, even if it is unimodal. The scenario gives a simple 
and cohesive explanation for the key observations, including (a) the overall shape of color 
histograms, (6) the number ratio of blue and red GCs as a function of host galaxy luminos- 
ity, and (c) the peak colors of both blue and red GCs as a function of host luminosity. If 
the bona-fide shape of the color-metallicity relationship is highly inflected, what has been 
thought to be the MDFs of GC systems may deviate significantly from the true distributions. 
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In this paper, we present an alternative way of resolving the long-standing discrepancy in 
the MDFs between GCs and halo stars in bright elliptical galaxies. The nonlinear conversion 
from metallicities to colors (Paper I) should not be irreversible, and here we try to inverse- 
transform color distributions of GCs into metallicity distributions using the nonlinear OMRs. 
Section 2 presents the pros and cons of the nonlinearity of OMRs, on which the present work is 
based. Section 3 applies the nonlinear color-to-metallicity conversion to the actual GC color 
distributions, and examines the inferred [Fe/H] distributions for M87 and M84 (^^ 3.1) and 
for the 100 early-type galaxies in the ACS Virgo Cluster Survey (ACSVCS, ICote et al.ll2004j ) 
(§§3.2). Section 4 compares the inferred GC MDFs with the MDFs of spatially-resolved 
constituent stars of nearby galaxies (§§4.1 and §§4.2) and the MDFs from a simple chemical 
enrichment model of galaxies (§§4.3). Section 5 discusses the implications of our results on 
the color bimodality issues (§§5.1), addresses the question of whether the GC formation is 
coupled with the bulk formation of the stellar population of host galaxies (§§ 5.2), and finally 
presents our view on the formation and evolution of GC systems and their parent galaxies 
(§§5.3). 



2. NONLINEARITY OF COLOR-METALLICITY RELATIONS: 

PROS AND CONS 

The core of this work is that the MDFs derived from GC optical color distributions are 
similar to those of constituent halo stars in galaxies. This work is based on the nonlinearity- 
CMR hypothesis (Paper I), that has been a target of dispute after its announcement. The 
issue is important enough that we devote a section to discuss it. 



2.1. Observed and Predicted Color-Metallicity Relations 



Simple linear conversion of photometric colors is frequently used for estimating metal- 
licities for large samples of extragalactic GCs. This is a reasonable first-order assump- 
tion for obtaining mean metallicities, but for investigating the detailed structure of the 
MDF, including possible subpopulations, the form of the CMR must be known to higher 
order. However, the best empirical color- metallicity calibrations currently available (e.g.. 



Peng et al 



2011 



2006 



Lee et al.ll2008l : iBeaslev et al.l 120081 : ISinnott et al.ll2010l : IWoodlev fc Harris 
Alves-Brito et al.l 1201 ll ) exhibit notable observational scatter. Moreover, compared 



to tens of thousands of GCs in a typical giant elliptical, the calibration samples are still 
relatively small and sparsely populated at the high-metallicity end. Larger samples of high- 
quality spectroscopic metallicities are needed to establish the precise forms of CMRs. Such 
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samples would implicitly include any correlations with age or other parameters, and would 
provide strong constraints on the theoretical models. 

In general, the slope of the dependence of a given photometric color on the logarith- 
mic metallicity [Fe/H] will change as a function of metallicity, and several recent stud- 
ies have found departures from linearity. For instance, it has been known for decades 
that the color of the gi ant branch in Galactic GCs is a nonlinear function of [Fe/H] (e.g.. 



Michel fc Smith 



Harris fc Harris 



9841) . iRichtlen (120061 ). using the observed C — Ti versus [Fe/H] relation of 



(120021 ). showed that the shape of the color distribution would differ from 



that of the M DF, and a r ion-p eaked metallicity distribution could result in a bimodal color 



distribution. iLee et al.l (120081 ) fitted nonlinear relations for C — Ti color versus [Fe/H], 



and found that the inferred MDF c hanged significantly depending on the adop t ed rel ation. 
Blakeslee. Cantiello. fc Pengl (|2010[ ) fitted a quartic relation to the iPeng et al.l ( 120061 ) data 
and showed that the empirical fit produces bimodal g — z colors from unimodal MDFs. 

As the observations have improved, so have the models, and these also tend to predict 
nonlinear OMRs (e.g., Lee et al. 2002; Paper I; Cantiello & Blakeslee 2007). This is partly, 
but not totally, due to improved modeling of the horizontal branch. Kissler-Patig et al. 
(1998a) showed that the Worthey (1994) models predict a nonlinear relation bet ween [Fe/H] 



and V — /, with a wavy form qualitatively similar to that found empirically by iPeng et al. 



(120061 ) . although the inflection occurs at higher metallicity because the colors of these models 
are generally too red at the high-metallicity end (see Blakeslee et al. 2001). This is interesting 
because the Worthey models do not realistically model the horizontal-branch morphology, 
but treat it as a red clump near the giant branch with a position that varies according to 
age. The Lee et al. (2002) model colors clearly showed nonlinear behavior as a function 
of metallicity even without any horizontal-branch component, but the nonlinearity in the 
optical colors was more pronounced with the horizontal branch. 

Despite the advances, more work is needed, especially on the behavior of the horizontal 
branch in extragalactic GC systems, as this is a complex multi-parameter problem. Metallic- 
ity is the primary factor governing horizontal-branch temperature, and the transition occurs 
in a nonlinear way at intermediate metallicities (e.g., Lee et al. 1994). However, variations 
in other parameters, including age, helium content, and central density, can create signifi- 
cant scatter in horizontal- branch morphology at a given metallicity (e.g., Sandage & Wildey 
1967; Zinn 1980; Stetson et al. 1996, Sarajedini et al. 1997; Buonanno et al. 1997; Sweigart 
& Catelan 1998; see also the recent discussions by Yoon et al. 2008; Gratton et al. 2010; 
Dotter et al. 2010). There are also intricate intercorrelations among these parameters, as well 
as correlations with GC mass. For instance, the color-magnitude relation found among the 
brightest GCs in external systems (e.g. Harris et al. 2006; Strader et al. 2006; Mieske et al. 
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2006, 2010; Peng et al. 2009) is likely the result of a mass-metallicity relation. In Galactic 
GCs, the presence of an extreme blue horizontal branch correlates strongly with GC mass 
(Lee et al. 2007), and is likely related to the presence of helium-enhanced subpopulations 
(e.g., Norris et al. 2004; Lee et al. 2005b; Piotto et al. 2005; Yoon et al. 2008; Han et al. 2009; 
Gratton et al. 2010). More massive GCs also have smaller half-light radii and higher central 
densities (van den Bergh 1996). Finally, there is some evidence for correlations between 
age and metallicity in both the Galactic and extragalactic GC systems (Puzia et al. 2005; 
Beasley et al. 2008; Marm-Pranch et al. 2009), but the degree of correlation must depend on 
the formation history of the GC system. 

To summarize, the best current data indicate that optical CMRs for GC systems are 
nonlinear, and multiple sets of models support this general result. However, exactly how 
colors vary with metallicity depends on ages, the age-metallicity relations, and variations 
in other parameters such as a-element and helium contents. Thus, a complete picture of 
the multi-band color-met allicity behavior in extragalactic GC systems would include an 
understanding of the interplay of the various stellar population parameters. Lacking more 
stringent observational constraints, we consider only some simple, reasonable assumptions 
on the standard stellar parameters, and for the present work, we use our best predicted 
CMRs. 

2.2. Assessing the Evidence for Metallicity Bimodality 

2.2.1. The Milky Way 

The GC system of the Milky Way Galaxy follows a bimodal MDF (Zinn 1985). The 
bimodality is confirmed by the radial number density profiles of metal-poor and metal-rich 
GCs and their orbital characteristics. The metaUicity distribution of GCs and field stars in 
our Galaxy are known much more accurately than those in any other galaxy. For instance, 
the metal-poor GCs in the Milky Way have formed with an efficiency ~ 20 times greater than 
the metal-rich GCs with respect to their associated stellar populations. Although these were 
well known, there was little expectation that it would be a general property of much larger 
GC systems in giant ellipticals. Such galaxies contain 10-100 times as many GCs as the 
Milky Way, and reside predominantly in cluster environments where they likely experienced 
very different evolutionary histories (e.g., Peng et al. 2008); thus, the analogy with our own 
Galaxy was unclear. 

In the context of the Toomrc (1977) idea that elliptical galaxies are the remnants of dissi- 
pationally merged spirals. Ashman & Zepf (1992) discussed bimodal MDFs with a population 
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of higher metaUicity GCs forming in the major merger. Nonetheless, given the simphcity 
of that model, and the advances in hierarchical structure formation theory, it came as a 
surprise when most giant ellipticals exhibited bimodal GC color distributions (interpreted 
as linearly reflecting metallicity), and other scenarios were proposed to account for the bi- 
modality (e.g., Forbes et al. 1997; Cote et al. 1998; Beasley et al. 2002; Kravtsov & Gnedin 
et al. 2005; Muratove & Gnedin 2010). However, the presumed bimodality of the MDF is 
taken as an input constraint in these models, rather than being a clear prediction (apart 
from perhaps the original Ashman & Zepf scenario, which is difficult to reconcile with the 
more complex assembly histories of ellipticals found in cosmological simulations). In light 
of this, it is worth reexamining the direct evidence for bimodal MDFs in elliptical galaxies, 
and other galaxies with comparably large GC systems. 



2.2.2. M31 

M31 is the largest galaxy in the Local Group and contains over 450 confirmed GCs, three 
times as many as the Milky Way, and hundreds of additional candidates (Galleti et al. 2004; 
Huxor et al. 2011; Caldwell et al. 2009, 2011). Barmby et al. (2000) studied the optical and 
near-IR color distributions for a large sample of M31 GCs. Unlike in most giant ellipticals, 
the colors did not appear bimodal, but Barmby et al. suggested that errors in reddening 
and photometry could "wash out" the bimodality. A KMM analysis of the V — K colors 
favored a double Gaussian model over a single Gaussian with 92% confidence (< 2cr), but 
although the distribution appeared asymmetric, it did not show two distinct components as 
in the Milky Way. The result was similar for the sample of ~ 160 GCs with spectroscopic 
metallicities compiled by these authors (see their Figure 19). 

Galleti et al. (2009) presented a homogeneous set of metallicities from Lick index mea- 
surements for 245 GCs. Again, they found that multiple Gaussian component models were 
favored because the MDF is broad and asymmetric, but it lacks the two distinct metallicity 
peaks seen in the Milky Way (see their Figure 15). when the sample was restricted to M31 
GCs with errors < 0.3 dex, there was no appearance of bimodality. Galleti et al. (2009) con- 
clude, "The MD[F] of M31 GCs does not present any obvious structure like the bimodality 
encountered in the GC systems of the MW [Milky Way]. Nevertheless, the distribution for 
M31 clusters does not seem to be well represented by a single Gaussian distribution.... While 
clearly not conclusive, the above analysis suggests that there may be actual structures in the 
MD[F] of M31 GCs." 



In the most recent study of the M31 GC MDF, ICaldwell et al.l ( 120111 ) present high 



signal-to- noise spectroscopic data on the M31 GC system. Metallicities were estimated using 
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a calibration of Lick indices with [Fe/H] provided by Galactic GCs. Although Caldwell et al. 
sample does not include many outer-halo metal-poor GCs that would increase the significance 
of the metal-poor side of the MDF, the metallicity distribution of over 300 old GCs has 
a significant population of intermediate-metallicity GCs and is not generally bimodal, in 
strong distinction with the bimodal Galactic GC distribution. The MDF shows a broad 
peak, centered at [Fe/H] = —1, possibly with minor peaks at [Fe/H] = —1.4, —0.7, and 
—0.2, suggesting that the GC systems of M31 and the Milky Way had different formation 
histories. Given the complex accretion history of M31 (McConnachie et al. 2009), it is not 
surprising that the M31 GC MDF would possess significant structure, but the best current 
data do not present evidence for bimodality in the M31 GC system. We note in passing that 
the Balmer absorption lines (H/3, H7 and B.6) theoretically have nearly the same response to 
horizontal-branch stars as optical broadband colors in their index-metallicity relations (e.g., 
Lee et al. 2000 ; Chu l et al. 2011, in prep.; S. Kim et al. 2011 (Paper IV)). Remarkably, 



Caldwell et al.l ( l201ll )'s sample shows strong nonlinearity in the Balmer lines vs. metal 
line (<Fe>) relations (see their Figure 10), and, as a result, exhibits clear Balmer strength 
bimodality (see their Figure 6) which is a close analogy with optical color bimodality. 



2.2.3. Cen A and the Sombrero 

The case for MDF bimodality is better in NGC5128 (Cen A), an SO pec galax}0 at 
the center of its own small group (Karachentsev et al. 2007). Beasley et al. (2008) present 
spectroscopic metallicities for 207 GCs in this galaxy. The resulting MDF is skewed towards 
high metallicities and apparently has three closely spaced peaks (see their Figure 5), in 
contrast to the two well-separated peaks in the Milky Way. This difference in MDF structure 
likely reflects the very different accretion histories. A very similar MDF was found by 
Woodley et al. (2010) in a spectroscopic study of 72 NGC 5128 GCs. In this case, a unimodal 
distribution provided statistically the best flt, but their [MgFe] index distribution was better 
fltted with a double Gaussian model. In both studies, the NGC 5128 GC metallicities and 
[MgFe] values lack the sharply bimodal appearance of the optical colors, especially of B — V 
andV-I (Peng et al. 2004b). However, Spitler et al. (2008) flnd that the optical- [3. 6/im] IR 
color distributions for 146 NGC 5128 GCs are distinctly bimodal, providing good evidence for 
MDF bimodality. They also found that similar data for a smaller sample of GCs in NGC 4594 
(the Sombrero) did not provide signiflcan t evid ence for or against MDF bimodality in that 



galaxy. More recently, lAlves-Brito et al.l (120111 ) present a spectroscopic MDF for over 200 



GCs in this galaxy, which is bimodal with peaks at [Fe/H] ~ —1.4 and —0.6. 



1 u 



NED homogenized morphology," |http: / /nedwww. ipac.caltech.edu/ 
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2.2.Ji-. Giant elliptical galaxies 

Various studies of GC metallicities in giant ellipticals from spectroscopic and near- 
IR/optical photometric data (Puzia et al. 2002, 2005; Cohen et al. 2003; Strader et al. 2007; 
Hempel et al. 2007; Kundu & Zepf 2007) are discussed in detail by Blakeslee et al. (2010). 
We summarize here and refer the reader to that work for the full discussion. 

Despite the groundbreaking nature of many of these studies, the spectroscopic results 
tend to be limited by the sample sizes (< 1% of the population), sample definitions, or 
sensitivity to the treatment of the data. Thus, Blakeslee et al. (2010) conclude that the direct 
evidence for metallicity bimodality in giant ellipticals is weak. For example, the bimodality 
reported by Strader et al. (2007) in the sample of 47 M49 GC metallicities derived from 
the Cohen et al. (2003) Lick index measurements depended on the calibration above solar 
metallicity. Despite the very pronounced color bimodality, there was no significant evidence 
for bimodality in the metallicities reported by Cohen et al. (2003), which extended to higher 
metallicities, rather than forming a clump near the solar value as in Strader et al. (2007). 
Given the uncertainty in the high- metallicity calibration (see discussion in Cohen et al.), 
and the sample limitations (0.6% of the GC population, observed with two spectroscopic 
masks), the issue remains unresolved for M49. More recently, Foster et al. (2010) have 
studied the Call triplet (CaT) feature in a sample of 144 GCs in the Eridanus giant elliptical 
NGC 1407, which has prominent optical color bimodality. In Galactic GCs, the CaT index 
is linearly related to metallicity, at least for [Fe/H] < —0.4 (Armandroff & Zinn 1988), 
and models indicate that it is very insensitive to age (Vazdekis et al. 2003). Foster et al. 
(2010) find that bright GCs near the peaks of the color distribution have very similar CaT 
strengths, indicating very similar metallicities despite the wide separation in color space. 
The distribution of CaT-derived metallicities does not appear bimodal, but because of its 
asymmetry, it is better fitted by a double Gaussian model at the 2a level. However, if this 
result were interpreted as MDF bimodality, then the components would differ significantly 
in amplitude, width, and position from those implied by the optical colors; i.e., it would be 
a different bimodality. Foster et al. remark that if the metallicities are taken at face value, 
then the very different color and metallicity distributions could be reconciled by a nonlinear 
CMR causing a unimodal MDF to appear bimodal in color space. 



In a recent near-IR/optical photometric study, iKundu fc Zepj (120071 ) p resented the I—H 



color distribu tion of 80 GCs in M87, which shows bimodality. More recently, IChies- Santos et al 



teOlOl . l2011al Jbl) present (optical - near-IR) colors for the GC systems in 14 early-type galax- 



ies, and find that the bimodalit y becomes less evident in g — Kg if compared to g — z and 



even less pronounced in z — Kg. IChies-Santos et al.l ( 120101 ) point out that the disappearance 



of bimodality in these colors while evident in the optical g — z color could be attributed to 
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a nonlinear-CMR effect altfiougfi tfie observational uncertainties could also account for it. 
Finally, we note that our new studies combining ACS and WFC3/IR data for NGC 1399 in 
the Fornax galaxy cluster (Blakeslee et al. 2011, in prep.) and Subaru/MOIRCS near-IR and 
CTIO optical data for M60 and NGC 4365 in the Virgo (S. Kim et al. 2011, in prep.) find 
independently that the I — H and I — Ks color distributions are not significantly bimodal, 
despite the strong bimodality in the optical colors of the same sample. 

Overall, the available data on the MDFs of GCs in giant ellipticals are at best ambigu- 
ous. In cases where spectroscopy - the more direct measures of metallicity than colors - 
suggest bimodality, it is much less apparent than the dramatic double-peaked histograms of 
colors for the same galaxies (e.g., Peng et al. 2006). Thus, at least some of the observed 
color bimodahty is likely due to nonlinear behavior of colors with metaUicity. As discussed 
above, there is empirical evidence for such nonlinearity having a form that tends to pro- 
duce bimodal color distributions. Pipino et al. (2007), considering the results of Puzia et al. 
(2005), also concluded that color- metallicity nonlinearity would help significantly in recon- 
ciling the photometric and spectroscopic data. The question of the relative importance of 
color-metallicity nonlinearity and metallicity bimodality in producing the observed GC color 
distributions remains open. In the meantime, besides the necessity of better spectroscopic 
samples, it is worthwhile to explore the possibility of deriving MDFs from the observed col- 
ors under the assumption of nonhnear CMRs. The following sections use the latest stellar 
population models to invert the colors for very large photometric samples and to examine 
the implications of the resulting MDFs. 



3. COLOR AND METALLICITY DISTRIBUTIONS OF 
EXTRAGALACTIC GLOBULAR CLUSTER SYSTEMS 

Our main objective is to investigate the MDFs for GC systems when inferred from the 
nonhnear CMRs. Compared to the metallicity-to-color conversion shown in Paper I, the 
inverse-conversion from colors to metallicity is more susceptible to the inevitable incom- 
pleteness of current population synthesis models. With a theoretical CMR that is somewhat 
incorrect in the color direction, for example, the metallicity-to-color conversion will still give 
color distributions with the correct shape, but the inverse- conversion will yield erroneous 
metallicity distributions. Moreover, the inverse-conversion from colors to metallicity may 
be hampered by the varying observational uncertainties depending on the colors of inter- 
est. With these caveats in mind, however, careful inverse-conversions may shed light on 
the structure of the GC MDFs, including possible subpopulations. In this section, we ap- 
ply the transformations of Paper I to the color distributions of GC systems and present 
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their inferred MDFs for M87 and M84 (§§3.1) and for 100 early-type galaxies imaged in the 
ACSVCS (§§3.2). 

3.1. Globular Cluster Systems in M87 and M84 

In this section, we present the results of our multiband photometry for GC systems of 
M87 (NGC 4486) and M84 (NGC 4374). We have selected M87 and M84, giant eUiptical 
galaxies in the Virgo cluster, because they both have GC systems with confirmed color 
bimodality in g — z and are among very few elliptical galaxies with deep w-band observations 
available. We refer the reader to Yoon et al. (2011, hereafter Paper II) for greater details 
on the multiband photometry of the M87 GC system in the context of nonlinear CMRs. 
The archival F336W images of HS T/WFPC2 and HS T/WFC3 were used to obtain tiFsseiy 
for GC candidates in M87 and M84, respectively. Our upssew-^Sind catalogs were matched 
with ACS/WFC gpAibw- and zpsbOLP-^duud photometry of Jordan et al (2009). We hereafter 
refer to up^sQw, QPAihWi and zf%5qlp mags as u, and z, respectively. Jordan et al (2009) 
selected bona-fide GCs with their magnitudes, g — z colors, and sizes. We further employed 
color cuts in the w-band colors to filter out contaminating sources, especially background 
star-forming galaxies. We used 591 GCs in M87 and 306 GCs in M84 that have reliable u, 
g, and z measurements in common. The samples are w-band limited. 

The merit of the multiband observations is clear: Since the form of CMRs hinges on 
which color is used, the shape of the color distributions varies significantly depending on 
the colors in use. Hence a comparative analysis of the GC MDFs that are independently 
obtained from distributions of different colors will put the nonlinearity hypothesis to the 
test, as proposed in Paper 1 and Paper 11. Among other optical colors, the -u-band related 
colors (e.g., u — g and u — z) are theoretically predicted to exhibit the most distinctive 
CMRs from other preferred CMRs (e.g., for g — z), and thus the most adequate to the 
task. Furthermore, the li-band colors are significantly less affected by the variation in the 
horizontal-branch mean temperature, having less infiected, "smooth" CMRs than g — z lor 
given ages. Therefore, for instance, the conversion from u — g color distributions to MDFs 
via the {u — g)-\Fe/Y{\ relation should be more straightforward than the case of g — z. The 
reason why the CMRs for ti-band colors are less inflected than the g — z CMR is two-folded: 
(a) the integrated -u-band colors of main-sequence and rcd-giant-branch stars are smoother 
functions of metallicity compared to g — z, and (6) the ?i-band colors arc less sensitive to the 
horizontal-branch temperature variation, which is due to the fact that the blueing effect of 
the optical spectra with increasing horizontal-branch temperature is held back by the Balmer 
discontinuity where the li-band is located ( Yi et al. 2004) . Such properties make the w-band 
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colors good metallicity indicators for a wide range of age, and the u-hand color distributions 
are expected to be significantly different from distributions of other optical colors such as 
g — z, V — I, and C — Ti. See Paper II for detailed discussion on the w-band colors as a tool 
to probe the nonlinearity of CMRs. 

Figure 1 shows the observed color distributions of the M87 GC system from u, g, and 
z photometry and the inferred MDFs. First, Figures la and Id present the g — z vs. u 
and u — g vs. u diagrams, respectively. Figures 16 and le present the g — z and u — g 
color distributions, respectively. The g — z distribution of the M87 GCs unambiguously 
displays two peaks around g — z = 1.0 and 1.4. In contrast, the u — g distribution for the 
identical sample does not appear to have clear bimodality. One may argue that the larger 
observational uncertainties in tx-band weaken the bimodality. Table 3 shows that the typical 
photometric error oi u — g is 2.2 times larger than that oi g — z for the entire sample, but 
at the same time the ranges spanned by the colors are /S.{g — z) = 1.1 mag and A{u — g) = 
2.1 mag, that is, the baseline of w — is 1.9 times longer than that of g — z. As a result, the 
relative sizes of error bars are {g — z : u — g) = (1.0 : 1.2). In a relative sense, the errors 
in the two colors are quite comparable to each other. Moreover, Paper II shows that the 
u — z color, which has smallest relative errors, still exhibits weaker bimodality in the color 
distribution compared to the g — z distribution. It is, therefore, not likely that bimodality 
in the u — g distribution of M87 GCs is simply blurred by larger observational errors in the 
M-band. 

The variation in the histogram shape for different colors may suggest that the form of 
the CMRs varies depending significantly on the colors in use. This is shown in Figures Ic and 
1/, along with our predictions from the Yonsei Evolutionary Population Synthesis (YEPS) 
models (Chung et al. 2011; Yoon et al. 2011, in prep.). In Figure Ic, the g — z colors are 
shown as a function of [Fe/H] for GCs in the Milky Way (open circles), and M49 and M87 



^The models in this study are constructed using the Yonsei Evolutionary Population Synthesis (YEPS) 
code. The YEPS model gen e rates (a) synth e tic col o r-magnitude diagrams for individual stars (see , e.g 



Lee. Demaraue. fc Zinnlll994 llee et al.|[l999l boOSbl: IRcv et all [ioOli IYoou fc Ledbooj IYoou et al 



200 



Han et al.lboOOl) and (b) synth etic integrated spectra for colors and a.bsorption indices of simple and coniposite 



stellar populations (see, e.g., Lee et al.j 2005a : Park fc Lee 1997 : Lee. Yoon. fc Lee 2000 : Rev et al 



20071 l2009l:E<:avirai et al.lbooaboOTaiffilRee et al.ll2007l: IYoou. Yi. fc Leell200dl2009HYoon fc Chundboogi 



2005 



Spitler, Forbes, fc Beaselvl 120081 iMieske et all 120081 IChoi. Goto, fc Yoonll2009l : ICho et al.ll201ll: lYoon et al 
20111 ) ■ One of the main assets of our model is the consideration of the systematic variation in the mean color 
of horizontal-branch stars as functions of metallicity, age, and abundance mixture of stellar populations. 
The standard YEPS model employs the Yonsei- Yale stellar evolution models (Y. Kim et al. 2002; Han et al. 
2011, in prep.) and the BaSeL flux library (Westera et al. 2002). The spectro-photometric model data of 
the entire parameter space are available at http: / / web.yonsei.ac.kr /cosmic/data/YEPS.htm , 
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(filled circles and triangles). The references to the observed data used in the relations are 
summarized in Table 1. The fifth-order polynomial fit to our model data for 13.9-Gyr GCs 
is overlaid (thick solid line, sec Table 2). For a comparison, the straight grey line represents 
the linear least-squares fit to the data. Figure 1/ is the same as Figure Ic, but for u — g 
color. Open circles, blue and red filled squares represent GCs in the Milky Way, M87, and 
NGC 5128, respectively. The m - ^ colors of the GCs in the Milky Way and NGC 5128 
were obtained from their U — B colors via the equation, [u — g) — 1.014 [U — B) -\- 1.372, 
derived from model data for synthetic GCs with combinations of age (10 ~ 15 Gyr of 1 
Gyr intervals) and [Fe/H] (—2.5 ~ 0.5 dex of 0.1 dex intervals). U- and S-passband are 
relatively close to u- and ^-passband, respectively, and thus U — B has responses to the 
horizontal-branch morphology in a way that is very similar to u — g. Therefore, U — B are a 
good proxy to u — g, and the U — B vs. u — g relationship is best described by a linear fit over 
a range of ages and metallicities (Table 1). Guided by the model, the form of the observed 
CMRs appears to vary from the {g — 2;)- [Fe/H] relation to (u — 51)- [Fe/H]. One may argue, 
however, that current data appear to be fit both by the theoretical nonlinear relations and 
the empirical straight relations, as there is only a weak indication that the modeled relations 
are actually better fits to the g — z ys. [Fe/H] and u — g vs. [Fe/H] data. Given the current 
level of observational accuracy and inhomogeneity of the data, the purpose of our study is 
not to determine the exact shape of color-metallicity relationships, but to investigate the 
consequences and implications of the possible nonlinear transformations. 

We now consider the inferred GC MDFs. Figures Ig and Ih show the MDFs for the M87 
GC system. On the one hand. Figure Ig presents the GC MDFs converted from g — z (red 
histogram) and u — g (blue histogram) colors that are based on the traditional linear color- 
to-metallicity conversion (thin grey lines in Figures Ic and 1/), and thus are just replicas of 
their color histograms shown in Figures lb and le. Note that both the overall shape and the 
peak positions do not appear to agree between the GC MDFs from the two colors. On the 
other hand. Figure Ih presents the GC MDFs converted from g — z and u — g colors that 
are based on the improved inflected relationship between color and metallicity (thick black 
lines in Figures Ih and le). In contrast to Figure Ig, the inferred GC MDFs in Figure Ih 
are modified drastically to have a strong metal-rich peak with a metal-poor tail. The two 
histograms in Figure Ih are more consistent with each other in terms of their overall shape 
and peak positions than those shown in Figure Ig. We note that our stellar population 
models show that, for given input parameters, the absolute quantities of output are rather 
subject to the choice of model ingredients such as stellar evolutionary tracts and model flux 
libraries; the different choices can result in up to ~ 0.2 mag g — z and u — g variation among 
models and the inferred [Fe/H] values accordingly (up to ~ 0.5 dex). Hence, one should put 
more weight on the relative values of inferred GC MDFs, i.e., the overall morphology of the 
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MDFs and their unimodality. We, however, wish to emphasize that the typical GC MDF 
shape is obtained invariably from different colors, i.e., g — z and u — g for M87 GCs. 

Figure 2 is the same as Figure 1, but for the GC system in M84. First, Figures 2a 
and 2d present the g — z vs. u and u — g vs. u diagrams for the M84 GCs. Figures 2h 
and 2e present the g — z and u — g color distributions, respectively. As for the possible role 
of observational uncertainties in weakening bimodality of the u — g color distribution, the 
typical photometric error of w — is 1.7 times larger than that of g — z for the entire sample 
(Table 3), but at the same time the ranges spanned by the colors are A{g — z) = 1.1 mag 
and A{u — g) = 1.8 mag, that is, the baseline of — is 1.6 times longer than that of g — z. 
As a result, the relative sizes of error bars are [g — z : u — g) = (1.0 : 1.1). In a relative sense, 
the errors in the two colors are comparable. It is, therefore, not likely that bimodality in the 
u — g distribution of M84 GCs is simply blurred by larger observational errors in the u-hand. 
Figures 2c and 2/ are the same as Figures Ic and 1/, respectively. Finally, Figures 2g and 
2h show the inferred MDFs for the GC systems in M84. Again, the inferred GC MDFs 
in Figure 2h have a strong metal-rich peak with a metal-poor tail. The two histograms in 
Figure 2h show a better agreement with each other in terms of their overall shape and peak 
positions than those in Figure 2g. 

In this section, we have obtained multi-band colors of GCs in the two representative gi- 
ant elliptical galaxies, M87 and M84, and examined their color and metallicity distributions. 
We have found that the distributions of different colors can be transformed into unimodal 
metallicity distributions that are strongly peaked with a broad metal-poor tail. The im- 
plications of the typical shape of the inferred GC MDFs and its similarity to those from 
chemical evolution models and field-star observations (see §4) will be discussed in §5. We 
note, however, that the similarity itself between the GC MDFs from multiband colors does 
not necessarily represent evidence that the model is correct. Whether the similar MDFs 
from various colors can be taken as evidence for the nonlinear-CMR scenario for the color 
bimodality is a sufficiently involved issue and fully explored in Paper II. 



3.2. Globular Cluster Systems in the ACS Virgo Cluster Survey 

Motivated by the findings above in §§3.1 for the individual galaxies and to avoid possi- 



(Cote et al. 


2004; 


Pens; et al. 


2006; 


Jordan et al. 


2009) 



transformation scheme to ~10,000 GCs in ACSVCS, the largest and most homogeneous pho- 
tometric database of extragalactic GCs currently available. Figure 3 presents the observed 
color distributions and inferred MDFs of GC systems in bins of host galaxy luminosity. In 
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Figure 3a, we show the observed color histograms of G C systems for seven bins of host galaxy 



magnitude. The data are the same as in Figure 6 of iPeng et al.l (120061 ) and here we listed 
the data in Table 4. The histograms are normalized by the GC number at their blue peaks, 
and multiplied by constants, C, for clarity. The magnitude bins are 1 mag wide and extend 
from Mb ~ -21.5 (-22 < Mb < -21, red, C = 1.0) to ~ -15.5 (-16 < Mb < -15, purple, 
C = 0.4). A Gaussian kernel of a{g — z) = 0.05 is applied. Clearly, the histograms appear 
bimodal or asymmetric across the entire luminosity range, with all GC systems containing 
blue peaks and with more prominent red peaks in brighter hosts. 

A close scrutiny of Figure 3a reveals the tendency of the dip positions (and blue peak 
positions) in the color histograms to become progressively bluer as the host luminosity de- 
creases. In the context of nonlinear CMRs, this can be explained if GCs in fainter galaxies are 
slightly younger than those in brighter galaxies. This is because at younger ages, the blue 



horizontal branch develo ps at lower metallicity ( iLee. Demarque. fc Zinrull994j : lYoon et al. 



20081 : iDotter et al.l 120101 ) and, as a consequence, the predicted colors of the quasi-inflection 
points along the CMR move systematically towards the blue at younger ages. This ef- 
fect is demonstrated in Figure 3b by an example set of the YEPS model predictions with 
At (brightest-faintest) = 3 Gyr. In this example, ages range from 10.5 Gyr (purple) to 13.5 
Gyr (red) by equal intervals of 0.5 Gyr. The YEPS g — z data for 9 Gyr to 14 Gyr by steps of 
0.5 Gyr are given in Table 2. The thin straight dotted line represents the linear least-squares 
fit to the data shown in Figures Ic and 2c. 

Figures 3c- 3/ present the results of the four different color-to- metallicity transforma- 
tions. Like the color distributions in Figure 3a, the inferred GC MDFs are displayed for 
seven bins of Mb from Mb = —22 to —15 in steps of 1 mag. To obtain these MDFs, we 
converted the color of each GC into [Fe/H] using the {g — 2;)-[Fe/H] relations from the YEPS 
model. Each panel makes a different assumption on the systematic age sequence from the 
faintest host bin to the brightest. The modeled ages for host luminosity bins are shown in 
the insets of Figures 3c- 3/, with age differences. At = 0, 1, 2, and 3 Gyr, respectively, 
between the faintest {Mb — —15.5, purple) and brightest (—21.5, red) bins. The brightest, 
oldest [Mb — —21.5, red) bin is set to be 13.5 Gyr. The color distribution shown in Figure 
3a was transformed to MDFs via the model CMRs of the corresponding ages. 

Figure 3c presents the case in which the age is assume to be constant at 13.5 Gyr re- 
gardless of host luminosities between Mb — —15.5 (purple) and —21.5 (red). The GC MDFs 
for the luminous {Mb < —17) host bins (the first to fifth brightest bins) in particular are 
strongly peaked with a broad metal-poor tail. As Paper I suggested, the strong bimodality 
seen in the GC color distribution of luminous galaxies is not evident in the MDF once trans- 
formed by their wavy CMR. The MDFs in the faintest two bins, however, have broad peaks. 
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This is likely because the less luminous galaxies primarily have blue GCs and their colors 
are largely bluer than the inflection point in the OMR under the constant age assumption, 
i.e., At(brightest - faintest) = Gyr. Paper I mentioned that the position of the inflection 
is bluer for younger stellar populations and so if fainter galaxies host GCs younger than 
brighter counterparts, then a different CMR may apply. However, in this naive use of the 
single CMR, not all inferred [Fe/H] distributions appear to be unimodal. 

Comparative analysis shows that the non-zero age difference of At(brightest - faintest) 
= 1 ~3 Gyr (Figures 3d-3f) results in the GC MDFs for all the host luminosity bins that flt 
better with skewed Gaussian distributions with metal-poor tails. The grey dotted histogram 
in each panel represents the inferred MDF for the brightest (i.e., Mb — —21.5) bin based on 
the simple straight fit in Figure 36, and is intended for comparison to the red sohd MDF. 
Obviously, the strong bimodality seen in the color distributions of luminous galaxies is no 
longer evident in the MDFs, once transformed via the inflected color-metaUicity relationship. 
The inferred GC MDFs for all seven host luminosity bins have the same characteristic shape 
- being sharply peaked with a broad metal-poor tail - across three orders of magnitude in 
the host galaxy mass. In addition, the mean [Fe/H] and peak position of the GC MDFs are 
a strong function of the host luminosity, in the sense that, for brighter host galaxies, the 
mean [Fe/H] increases and the peak gets redder. 

Compared to their relative ages, the absolute ages of GC systems are still less certain. 
To test the robustness of the result shown in Figure 3 against different absolute ages. Figure 
4 makes differing assumptions on the age sequence from the faintest host galaxies to the 
brightest. In this case, the center bin, i.e., the fourth brightest {Mb — —18.5, green) bin, 
is set to be 13 Gyr. The modeled ages for the host luminosity bins are shown in the insets 
of Figures 4c-4/, with age differences of At — 0, 1, 2, and 3 Gyr, respectively, between 
the faintest {Mb ^ —15.5, purple) and brightest (—21.5, red) bins. In Figure 46, another 
example set of our model prediction with At — 3 Gyr is shown. In this example, ages range 
from 11.5 Gyr (purple) to 14.5 Gyr (red) in equal intervals of 0.5 Gyr (Table 2). Figures 
4c -4/ show that the strong bimodality is no longer evident in the inferred MDFs, and the 
typical shape of the GC MDFs docs not depend upon differing age assignment from the 
faintest host galaxies to the brightest. It also holds true that the peak positions of the GC 
MDFs are a strong function of the host luminosity in that, for brighter host galaxies, the 
mean [Fe/H] increases and the peak gets redder. We also tested the stabihty of the results 
against the putative age dispersion in the examined data. The typical form of the GC MDFs 
persists with up to at — 2.5 Gyr, where the age spread among GCs in a single host luminosity 
bin is parameterized by a Gaussian dispersion, at- The data for the inferred GC MDFs in 
Figures 3d and 4d (they are identical) are given in Table 5. 
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If one could assume that the intrinsic shapes of the GC MDFs are similar in all host 
luminosity bins, our results might give a tantalizing hint of varying ages among GC systems 
in that GCs in more luminous parents are older than GCs in fainter galaxies. However, 
there is yet no observational support that the assumption is valid. Moreover, the true shape 
of CMRs suggested in this study is still unproven. Hence, the purpose of our simulations 
should be to determine neither the absolute age of the GC systems nor their exact pattern 
of age sequence from the faintest bin to brightest. Nevertheless, the possible age seque nce 
appears not inconsistent with the "galaxy downsizing" picture (e.g.. ICowie et al.lll996l ) in 
which brighter galaxies are observed to form earlier. If confirmed, this can be regarded as 
another indication of the common characteristics shared by stellar populations of GCs and 
halo field stars (see §4). 



4. COMPARISON OF GC MDF'S TO THOSE OF HALO FIELD STARS 
AND GALAXY CHEMICAL EVOLUTION MODELS 

We find that the bimodal GC color distributions commonly observed in luminous early- 
type galaxies are transformed into unimodal metallicity distributions that are strongly peaked 
with a broad metal-poor tail. Our sample includes GC systems in M87 and M84 (§§3.1) and 
in the ACSVCS galaxies (§§3.2). A key to comprehending the connection between GCs and 
halo field stars in galaxies is a direct comparison of their MDFs. This section compares the 
inferred GC MDFs to those of resolved field stars in nearby early-type galaxies (§§4.1 and 
§§4.2). Also, the GC MDFs are compared to MDFs produced by chemical evolution models 
of galaxies (§§4.3). 



4.1. Comparison of ACSVCS GCs to Spatially-resolved Halo Field Stars in 

Nearby Elliptical Galaxies 



The necessity of obtaining photometry of spatially-resolved field stars limits us to nearby 
galaxies. There are several nearby, relatively massive elliptical galaxies whose stellar MDFs 
have been measured, and our inferred GC MDFs can directly be compared to such stellar 
MDFs. Figure 5 gives a compariso n of the MDFs of the ACSVCS GCs to tho s e of resolved 



2007a b; Bird et al 



stars of nearby elliptical galaxies ( iHarris fc Harris! 120021 : iRejkuba et al.l l2005l : [Harris et al. 



20101 ). The stellar MDFs of individual galaxies were obtained from color- 



magnitude diagrams of red-giant stars whose colors are highly sensitive to their metallicities. 



In Figure 5, we plot the GC MDF shown in Figures 3(i and Ad (they are identical). The 
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data for the plot are presented in Table 5. The figure shows that the GC MDFs are similar to 
the MDFs of resolved constituent stars of nearby elliptical galaxies. In particular, the typical 
shape of the GC MDFs (colored solid line in each panel), characterized by a sharp peak with 
a metal-poor tail, is remarkably consistent with those of field stars in nearby galaxies. By 
contrast, GC MDFs obtained using the simple linear [Fe/H] vs. g — z relation (grey dotted 
lines) do not agree with the stellar MDFs. We note that the metallicity spread of GCs tends 
to be broader than that of stars. As indicated in Figures 3c- 3/ and Figures 4c-4/, the 
GC MDFs, on the whole, are broader than the stellar MDFs. This is likely due to the fact 
that each color histogram of the ACSVCS consists of GCs belonging to 3 ~ 20 diverse host 
galaxies, giving an ensemble character of each host luminosity bin. Moreover, the inferred 
MDFs become broader as the observational uncertainty in color is propagated to metallicity 
space. 

In addition to the noticeable similarities found in the shape of MDFs between GCs and 
stars, they share a common feature in that their mean metallicity increases with increasing 
host luminosity. As a result, the peak positions of GCs and stars are roughly coincident in 
each luminosity bin — at [Fe/H] ~ —1.0 (faint hosts) and —0.5 (bright hos ts). This, however 



should be taken with great caution because the rnean c olors of both G Cs (iDirsch et al. 



Jordan et al.l 



2004 



Tamura et al.ll2006l : iLee et al.ll2008l ) and field stars (IHarris fc Harris 



2003 



2002 



Rejkuba et al.ll2005l ) depend significantly on the sampled radial location in a galaxy. Their 
mean metallicities gradually decrease with projected radius, as seen in Figures 5d and 5e 
for the outer- and inner-halo stars of NGC 5128, respectively. Moreover, observ ations show 
that the mean GC color is bluer than stars within an elliptical galaxy as a whole (jPeng et al. 



20061 ). More importantl y, GCs are on average b 



at the same radii (e.g.. iForte, Strom, fc Strom 



Tamura et al. 2006: Lee et al. 



2008 



uer than the unresolved hght of the galaxies 



19811: 



Strom et al. 1981: Jordan et al. 2004 



Forte. Vega, fc Faifeiil2009l ). For old stellar populations 



such as those in elliptical galaxies and GCs, the observed bluer colors imply lower metallicities 
in all stellar models, thus leading to the conclusion that GCs are typically more metal-poor 
than field stars in a galaxy. Hence it may just be a coincidence that the GC MDFs line 
up with the field star MDF of the nearest elliptical galaxy. Despite these caveats, the peak 
and width of the MDFs of GCs and field stars are similar enough to suggest that both were 
formed in the same events, which built the major part of the galaxies. 

It should be addressed that the [Fe/H] histogram for NGC 3379 stars (Figure 5c), despite 
the overall resemblance to t he compared GC MD F, shows an excess of metal-poor stars with 
a bump at [Fe/H] ~ —1.2. [Harris et al.l (j2007bl ) found the metal-poor stellar halo in NGC 
3379 and provided an explanation for why their earlier studies only detected a metal-rich 
component in the same galaxy. The NGC 3379 HST field is at a distance of 12 effective radii, 
which is more than twice as far as the equally small fields of view for NGC 3377 (Figure 5b) 
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and NGC 5128 (Figures 5d and 5e). Studies of the stellar MDFs of large galaxies have been 
performed photometrically and concentrate on more central regions of the galaxy where the 
metal-rich field star population is concentrated. The radial bias seemed to under-sample 
metal-poor stars in the earlier studies. With the recent evidence of metal-poor halo stars 
emerging at great galactocentric distances, one should not assume that the stellar MDFs are 
unimodal at all positions. Therefore, this suggests that our interpretation should be more 
applicable to their main, inner parts of galaxies, less to the remote outskirts. 



4.2. Comparison between Spatially-resolved Halo Field Stars in Nearby 
Elliptical Galaxies and their Own GC systems 

As envisaged by the NGC 3379 case above, each galaxy has its own evolutionary history. 
So a more direct way to assess the similarity between GCs and stars is to compare between 
the MDFs of spatially-resolved field stars in a galaxy and its own GC system. There are four 
galaxies (M87, NGC 5128, NGC 3377, and NGC 3379) in Figure 5, for which both the GC 
MDF and the stellar MDF are currently available. Figure 6 gives the comparison for M87 
(the top row), NGC 5128 (the second row), NGC 3377 (the third row), and NGC 3379 (the 
bottom row). The stellar MDFs (grey histograms) are identical to those in Figures 5b, 5c, 
5d, and 5/. The left-hand panels present the inferred GC MDFs (empty histograms) based 
on the traditional linear color-to- metallicity transformations. By contrast, the right-hand 
panels show the inferred GC MDFs (empty histograms) based on the inflected relations from 
the YEPS model (Table 2). 



Firstly, for M87, we exploit g—z colors from the ACSVCS ( jPeng et al.ll2006l : I Jordan et al. 



20091 ) to derive the GC MDF. Since the inferred GC MDFs shown in Figure Ih are based on 
the M-band limited sub-sample, the ACSVCS g- and z-band data are more representative GC 
sample. Figure 6a shows that the MDF (empty histogram) for 1745 GCs obtained using the 
simple linear [Fe/H] vs . q — z relation ( shown in Figure Ic) exhibits a fundamental difference 
from the stellar MDF ( iBird et al.ll2010l ) measured in a similar region of the same galaxy. By 
contrast, in Figure 6e, the strong bimodality is no longer present once transformed by the 
inflected CMR (Table 2), and the GC MDF is very similar to that of the brightest host bin 
in the ACSVCS. As a consequence, the MDFs of GCs and field halo stars in M87 are similar 
in shape and line up remarkably well with each other. 



Secondly, for NGC 5128. |Peng. Ford, fc FreemanI fl2004al Jbh presented the CTIO Blanco 
4-m U-, B-, V-, and /-band photometry of the GC system. We use the B—I color distribution 
to derive the GC MDF because B — I is a. reasonable substitute for g — z. However, our 
result on NGC 5128 is not affected by the choice of colors. Figure 66 shows the MDF (empty 
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histogram) for 210 GCs that has strong bimodahty. On the other hand, in Figure 6/ we 
converted the B — I colors into [Fe/H] using the YEPS relation (Table 2). The strong 
bimodality seen in Figure 6fe is not evident in the MDF, and the strongly peaked GC MDF 
with a broad rnetal-p oor tail is in reasonably good agreement in shape with the stellar MDF 
jReikuba et aliboosh . 



An unavoidable uncertainty in this comparison is that the field stars are sampled from 
only one or two projected location(s), whereas the GCs cover the wider halo. For M87, the 
stars were sampled in an inner region (R ~ 10 Kpc) and the GCs were in the inner halo (< 
10 Kpc). For NGC 5128, the stars were sampled in an inner region (R ~ 20 and 30 Kpc) 
and the GCs were in the entire halo of the galaxy. This may partly explain the similarity 
of the peak positions between stellar and GC MDFs for M87, and the dissimilarity for NGC 
5128. Another warning for NGC 5128 is that, althou gh the majority of its GCs (~ 90 %) 
are known to be old (> 10 Gyr) (IBeasley et al.l 120081 ) . there should be a certain portion of 
young GCs for which one needs to apply a different color vs. metallicity relation to derive 
the GC MDF. Nonetheless, despite these sources of uncertainty, it is engrossing that the 
inferred GC MDFs based on the inflected color-to-metallicity transformations gives better 
matches with the fleld star MDFs of M87 and NGC 5128, compared to those based on the 
traditional linear relations. 



Thirdly, for NGC 3377, ICho et al.l mm presented the HST/ ACS g- and z-band pho- 
tometry of the GC system as part of a deep imaging study of 10 early-type galaxies in 
low-density environments. Figure 6c shows that the MDF for 157 GCs derived from the 
simple linear [Fe/H] vs. g — z relation has strong bimodality. On the other hand, in Figure 
6g we converted the g — z color into [Fe/H] using the YEPS relation. The strong bimodality 
seen in Figure 6c is not evident in the MDF, and the sharply peaked MDF with a broad 



metal- poor tail is in reasonably good agreement in shape with the stellar MDF (IHarris et al. 
2007313) • 



Lastly, for NGC 3379, IWhitlock. Forbes, fc BeaslevI (120031 ) and iRhode. fc Zepj (120041 ) 
carried out wi de-fleld photonietry o f GCs. They verifled earlier results that its GC population 
is quite small. iHarris et all J2007bh derived the GC MDF from Rhode & Zepf (2004) 's 5-/2 
histogram for 36 GCs with an improved empirical relation, [Fe/H] = 3.13 {B — R)o - 5.04. 
Figure 6d shows the GC MDF, which is highly concentrated toward the metal-poor side. If 
divided at [Fe/H] = —1.2 where the metal-poor bump of the filed-star MDF is located, the 
metal-poor GCs outnumber the metal-rich GCs by 25 to 11. As a result, for the metal-rich 
half of the MDF, the numbers of GCs are much too small. By contrast. Figure 6h displays 
the GC MDFs based on the YEPS relation (Table 2). With the GC MDFs of only 36 
clusters, it is not yet clear whether the MDFs of GCs and stars have the same shape. It is 
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interesting to note, however, that the GC metalhcities are now more evenly distributed, and 
the proportions of the two metalhcity subgroups, when divided at [Fe/H] = —1.2, became 
more comparable (21 : 15) than those in Figure 6d (25: 11). Interestingly, the inferred GC 
MDF in Figure 6h seems very similar to the NGC 3379 field-star MDF in the western half 
(farther from the galaxy center) of the ACS/WFC field (shaded histogram in Figure 14 of 
Harris et al. 2007b) 



4.3. Comparison to Chemical Emichment Models of Galaxies 



The new development in linking the observed GC colors to their intrinsic metalhcities 
leads to GC MDFs that are strikingly similar in shape to MDFs of resolved field stars in 
nearby elliptical galaxies. Both the inferred GC MDFs and halo stellar MDFs are charac- 
terized by a sharp peak with a metal-poor tail. In this section, we proceed to compare the 
inferred GC MDFs to chemical enrichment models of galaxies. 

In Figure 7 we compare the inferred GC MDFs for the second brightest {Mb — —20.5, 
orange) and the brightest {Mb — —2 1.5, red) galaxy luminos ity bins with the simple closed- 



box model of chemical evolution (e.g. JPagel fc Patchettlll975l ). We plotted the MDFs for this 



simple model using yields of [Fe/H] = —0.685 and —0.505, respectively. With the improved 
nonlinear relationship between color and metallicity, the general shape of GC MDFs is in 
remarkable agreement with that of galaxy chemical enrichment models. In contrast, the 
dotted line in each panel represents the inferred GC MDF for the corresponding bin based 
on the simple straight fit to the data. 

Despite their resemblance in the general shape, the inferred GC MDFs have fewer metal- 
poor GCs than the simpl e chemical model. This is akin to t he "G-dwarf problem" in the 
solar neighborhood (e.g., Ivan den Berghl Il962l : ISchmidtl Il963[ ). We note that the width of 



the chemical evolution model can be changed using various kinds of gas infall and stellar 
feedback, and an accreting-box model of chemical evolution yields narro wer distributions 
that p rovide a better match with the inferred GC MDFs. For instance, [Harris fc Harris 
fc002h remarked that the halo of NGC 5128 {Mb = -20.9) also suffers from the G-dwarf 
problem in that it lacks metal-poor stars compared to the simple model. They were able 
to fit the field star MDF using an accreting-box model of chemical evolution, producing a 
narrower distribution that is also a better match to the inferred GC MDFs in this study. 
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5. DISCUSSION 

There is substantial evidence that GCs are the remnants of star formation events in 
galaxies, and are linked to th e star formation, chemical enrichment, and merging histories 



of their parent galaxies (e.g., iMcLaughliru Il999l ) . However, ever since direct photometry of 



spatially-resolved constituent stars in a dozen nearby galaxies became possible thanks to the 
HST and large ground-based telescopes, the discrepancy between the MDFs of GCs and field 
stars has remained a conundrum. If GCs mirror field stars across the galaxy histories, the 
MDFs of GCs and field stars should be similar. Thus, the curious disagreement in metallicity 
has been interpreted in the context of highly decoupled formation and evolution histories 
between GCs and constituent stars of their parent galaxies. 

Current observational data and modeling point convincingly to nonlinear CMRs, which 
have significant implications for the interpretation of GC color distributions. We find that 
the strongly peaked [Fe/H] distributions inferred from nonlinear CMRs are qualitatively 
similar to the MDFs of field stars in the spheroidal component of nearby galaxies and to those 
produced by chemical evolution models of galaxies. However, whether the inferred GC MDFs 
represent the intrinsic, true ones is still unproven, and so it may be partly a coincidence. 
Nevertheless, if the MDFs obtained using stellar population models more closely represent 
the true GC MDFs, then this would change much of the current thought on the formation of 
GC systems and their host galaxies. The next two sections discuss what constraints would 
our findings pose on formation of GC systems (§§5.1) and their host galaxies (§§5.2). In 
Section 5.3, we present our view on the formation and evolution of GC systems and their 
parent galaxies. 



5.1. What Do the Inferred GC MDFs Imply? 

Remarkable progress has occurred over the past few decades in our understanding of 
extragalactic GC systems. One of the most important discoveries is that many galaxies 
show bimodality in their color distributions, leading to the notion that galaxies possess two 
distinct subpopulations of GCs. 

We showed, however, that the typical GC MDF shape derived from color distributions is 
unimodal and characterized by a sharp peak with a metal-poor tail. If confirmed, the inferred 
GC MDFs may appreciably reduce the demand for the separate formation mechanisms to 
explain the metal-poor and metal-rich division of GCs. We warn, though, that the sample 
used in this study is the GC systems in the Virgo galaxy cluster obtained from the HST 
ACS/WFC, WFPC2, and WFC3 observations. The field of view of ACS/WFC, for example. 
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covers galaxies' haloes within R ~ 0.6, 0.8, 7.2, 8.3, 8.9, 9.4, and 10.6 Re (in z-band, 
Ferrarese et al. 2006) for the seven host luminosity bins (from the brightest bin to the 
faintest) of the ACSVCS, respectively. Nevertheless, the dearth of metal-poor GCs in MDFs 
for inner spheroids of giant {Mb ^ —20) elliptical galaxies and nearly the entire spheroids 
of normal (—20 < Mb ^ —15) ellipticals suggest that the inheritance of metal-poor GCs 
via dissipationless accretion from dwarf satellites seems to be less significant than previously 
thought. 

On the outskirts of giant galaxies, accretion of metal-poor GCs from low-mass satellites 
and/or from surrounding regio ns may be an important channel for a galaxy to add GCs on the 



metal-poor part of GC MDFs flForte. Martinez, fc Muzziolll982l:IC6te. Marzke. fc Westlll998 



Cote. West, fc Marzkd l2002l : Masters fc AshmanI l2010l : iLee et al.l l2010af ). Recall, however 



that one of the main consequences of the nonlinear metallicity-color relations is that their 
steepness at the metal-poor end naturally creates a blue peak of GCs in color space, which 
is a direct cause of the conventional subpopulation of blue GCs. Therefore, even for the 
outskirts of giant galaxies in cluster environments, whether or not accretion of metal-poor 
GCs is solely responsible for blue peaks of color distributions is still an open question. The 
strongest evidence against accretion models is that the blue peak colors are correlated tightly 
to the host galaxy luminosity (Larsen et al. 2001; Strader, Brodie & Forbes 2004; Peng et al. 
2006). The metal-poor relation implies that metal-poor GCs, although they formed at very 
high redshift and were accreted later on, already "knew" which galaxy they would ultimately 
belong to, and thus weakens the accretion scenario for the color bimodality. Alternatively, 
the nonlinear metallicity-color relations scenario (Paper I) gives cohesive explanations for 
the observations that the mean colors of both blue and red GCs increase progressively for 
more luminous host galaxies. 

Further counter-evidence of the accretion model is the significant fraction of blue GCs in 
massive cluster galaxies in relatively lower-density regions (Peng et al. 2006) and massive field 
galaxies in isolation (Cho et al. 2011). In such environments, galaxies have few neighboring 
lower-mass galaxies, and it would be difficult to acquire many metal-poor GCs via accretion. 
Therefore, the accretion process seems more important when it meets the three conditions: 
(a) the outskirts (rather than the inner, main bodies) of (b) giant (rather than dwarf) 
ellipticals orbited by a large number of low-mass satellites in (c) cluster (rather than isolated) 
environments. In this regard, wide-field, multiband studies of GC systems in cluster and field 
environments are clearly needed. Wide-field photometry of nearby cluster galaxies in CTIO 
4-m U-hand (H. Kim et al. 2011, in prep.) and Subaru/MOIRCS NIR (S. Kim et al. 2011, 
in prep.), and a study of massive field galaxies in HST/ ACS g and z (Cho et al. 2011) are 
done or in progress to further investigate our alternative scenario. 
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Current hierarchical models of galaxy formation in the ACMD cosmology predict that 
several thousands of small building blocks were involved for the emergence of one single 
massive galaxy. This is significant because the extent of complexity may leave little room for 
the existence of just two GC subpopulations in each massive galaxy. Indeed, the unimodal, 
skewed MDFs arise naturally in an aggregate of a large number of protogalactic gas clouds 
from its virtually continuous chemical evolution through many successive rounds of star 
formation. The strongly peaked unimodal GC MDFs point to GC formation with a relatively 
short, quasi- monolithic timescale. Remarkably, the typical GC MDF shape emerges across 
three orders of magnitude in host galaxy mass. This suggests that the processes of GC 
formation and chemical enrichment are quite universal among a variety of GC systems. 



5.2. Do GC Systems Trace Star Formation in Galaxies? 



We also address the important issue of whether or not the formation of GCs is coupled 
with the bulk formation of the stellar population of host galaxies. Or, equivalently, does 
GC formation really mirror star formation in a galaxy? We have shown that the inferred 
GC MDFs agree reasonably well with the stellar MDFs of nearby galaxies and the MDFs 
produced by models of galaxy chemical evolution. The results suggest that the evolutionary 
histories of GC systems and their parent galaxies are strongly coupled, and thus share a 
more common origin and closer subsequent evolution than previously thought. 

An important aspect of the GC-host galaxy co-evolution issue concerns the GC-to-star 
offset, in the sense that GCs are on average more metal-poor than the unresolv e d light of an 
elliptical galaxy at the same radial location 



1981; Jordan et al. 2004; Tamura et al. 2006 



e.g., 



Forte. Strom, fc Strom 



1981 



Strom et al. 



Lee et al.lbood : iForte. Veea. fc Faifeni2009[ ). In 



light of the typical shape of GC MDFs proposed in this study, the lower value of the mean 
metallicity of GCs compared to that of field stars should not be attributed to the number 
excess of the metal-poor GCs. Alternatively it is more likely that, at a given radial location, 
the GC MDF on the whole is shifted toward the metal-poor side with a peak at a lower 
[Fe/H] value, compared to the field-star MDF. 

The GC-to-star offset would point to a picture in which GCs are the remnants of vigorous 
star burst events in the early stages of galaxy formation, and thus preferentially trace the 
major mode of star formation in galaxies. If so, GC formation was less prolonged than 
field star formation. Recent observati ons show that GCs were at least an order of magnitude 
more massive at birth than now (e.g.. lConroyll201ll ). and the large masses may have been the 



cause of the earlier truncation of their formation process than that of stars. Consequently, 
the chemical enrichment process of a GC system appears to have ceased somewhat earlier 
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than that of the field stellar population in each star formation episode. 



We further showed a possible age difference among GC systems, in that the GC systems 
in fainter galaxies are on average younger. We refer to this as "GC system downsizing." The 
GC system downsizing phenomenon appears to further support the similar nature shared 
by stellar populations of GCs and field stars. Interestingly, recent observations reveal that 
the GC systems in Milky Way satellites and some Milky Way GCs believed to be accreted 
from satellite dwarf galaxies are relatively younger than the majority of Galactic GCs (e.g., 
Marm-Franch et al. 2009) 

Galaxy downsizing is generally defined by more prolonged residual star formation in 
fainter galaxies. However, the "GC system downsizing" involves the idea that the first 
generation of GCs in fainter hosts are created later than the counterpart in brighter galaxies, 
and may provide further information on galaxy downsizing itself. Provided that GC systems 
became attached to the present host galaxies from the beginning, we can speculate that 
not only the first GCs but also first field stars were formed earlier in brighter galaxies than 
fainter galaxies. That is, massive galaxies today were likely where conditions first favored star 
formation, thus suggesting a prolonged epoch of galaxy formation in the universe. Combined 
with the fact that GCs and field stars in brighter galaxies have higher metallicities, a picture 
emerges in which the formation and accompanying metal enrichment of both GCs and halo 
stars seem to have started earlier and proceeded more rapidly and efficiently in massive 
galaxies, presumably in denser environments. 

It is also important to note that some galaxies have substructure in their stellar MDFs 
and are thus more complex than a smooth transition from a metal-r ich dominance to 



metal-poor with incr e asing radius. Examples include the Milky Wav (llbata et al. 



Maiewski e 


t al.| 


Koch et al. 


2008 



McConnachie et al. 



Maiewski et all 20031: lYannv et all l2003l: llvezic et all 120081) and M31 flKahrai et al. 



2001'; 
12006 : 



20091 



1). The stellar MDF of the e lliptical galaxy NGC 



3379 (Figure 5c) also shows a fine substructure ([Harris et al.l l2007bl ). The remote outer 
haloes of these galaxies are inhomogeneous in terms of surface density and show different 
metallicity distributions from the inner haloes. This indicates that the stellar populations in 
the outer haloes of galaxies ar e not well mixed, and in turn supports the build up of halo for- 



mation via satellite accretion (IForte. Martinez, fc Muzzidll982t ICote. Marzke. fc Westlll998 



Cote. West, fc Marzkel [200^ : [Masters fc Ashmanlboioh. which can have mixing time scales 
of a few Gyr in these outer regions (I Johnston. Spergel. fc Hernquistlll995l ). 



The evidence of the satellite accretion indicates that the outskirts of galaxies do not all 
build up their stellar populations in one single way. For instance, the NGC 3379 observation 
fits into a model in which its outskirts were formed by a combination of earlier dissipative 
mergers and later accretion of dwarf satellites. The remote halo of NGC 3379 has been 
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on the fringes of active, violent gas-rich mergers at the early epochs of galaxy formation, 
and later experienced dry accretion of stars and GCs from its metal-poor satellite dwarfs. 
Therefore, one should not assume that the stellar MDFs are unimodal at all positions, and 
our new interpretation would be more applicable to the inner, main bodies of galaxies than 
their remote outskirts. 



5.3. An Alternative View on the Formation and Evolution of GC Systems and 

their Parent Galaxies 

Our results may be an important step forward in resolving the long-standing disagree- 
ment between GCs and field stars, and in reconstructing the history of GC systems and their 
parent galaxies. Although the sample used in this study (the HST ACS /WFC, WFPC2, and 
WFC3 photometry for the GC systems in the Virgo galaxy cluster) confines our discussion 
to R < Re for giant ellipticals and < 10 Rg for normal ellipticals, our findings suggest that 
GC systems and their parent galaxies have shared a more common origin than previously 
thought, and hence greatly simplify theories of galaxy formation. The star formation and 
accompanying chemical evolution were virtually continuous via aggregates of a large number 
of protoclouds and via repeated gas-rich mergers with other galaxies, leading to the uni- 
modal, skewed MDFs of both stars and GCs. The metal enrichment of both stars and GCs 
proceeded more rapidly and efficiently in massive galaxies, resulting in more metal-rich stars 
and GCs in those galaxies. The observed radial metallicity gradients are understood if the 
chemical enrichment in the dense centers was more rapid and efficient than in the less-dense 
outskirts. The typical GC MDF shape emerges across three orders of magnitude in the host 
galaxy mass, suggesting that the processes of GC formation and chemical enrichment are 
quite universal among various GC systems, at lease for the inner, main spheroids of giant 
eUipticals and the nearly entire haloes of normal eUipticals. 

The histories of GCs and their host galaxies are reconstructed as follows. 

1. Giant ellipticals' inner, main haloes: The inner, main spheroids (R < 1 R^) 
of today's giant elliptical galaxies {Mb < —20) were first created in dense regions of 
the universe via dissipational mergers of a large number of protogalactic gas clouds 
(e.g., Searl & Zinn 1978). Their GC systems with the MDF peaks at [Fc/H] > -0.7 
(Figures 3 and 4) were formed together with stars in the galaxies. The GC systems 
from low to high metaUicities (i.e., both blue and red GCs), as a whole, do not need two 
separate mechanisms for formation. There is strong similarity between blue and red 
GCs in their mass function; for the less evolved high-mass part of the mass function, 
these are approximately power laws with indices of —1.8 to —2. Although power-law 
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distributions are a consequence of a variety of physical processes, their nearly identical 
power indices may support an identical formation history of blue and red GCs. The 
chemical enrichment of the galaxies was rapid, and the difference in the formation 
epochs of the metal-poor and metal-rich ends of MDFs should be small (< 1 Gyr), as 
evidenced by observations (e.g., Jordan et al. 2002). In this view, the metal-poor GCs 
in massive galaxies are the first generation of GCs in the universe. 

2. Normal ellipticals' entire haloes: The spheroids (out to i? < 10 -Re) of normal 
galaxies {Mb ^ —20) were created later than those of massive ones via self-collapses 
or mergers of protogalactic clouds that did not take part in massive galaxy formation 
early on. They were outside of tumultuous, dense regions, and have evolved rather 
independently from central massive galaxies. Their GC systems with the MDF peaks 
at [Fe/H] < —0.7 (Figures 3 and 4) were formed together with stars in the low-mass 
galaxies. As a result, the first generation of both stars and GCs in low- mass galaxies 
are younger than those of massive galaxies, indicating a prolonged epoch of galaxy 
formation in the universe. This picture is supported by recent observational evidence 
that the GC systems in Galactic satellites and some Galatic GCs believed to be accreted 
from satellite dwarf galaxies are relatively younger than the majority of Milky Way 
GCs (e.g., Marin- Franch et al. 2009). 

3. Remote outskirts of giant ellipticals in dense environments: The remote outer 
haloes of giant elliptical galaxies have been on the fringes of active, violent gas-rich 
mergers at the early epochs of galaxy formation, and later exp erienced dissipationless 



19821: 



Cote, Marzke. & West 


Lee et al. 


2010a 


) and/ 



Cote. West, fc Marzke 2002; Masters fc Ashman 



Bergond et al. 2007; Schuberth et al. 2008; Lee et al. 2010b; West et al. 2011). Such 
GCs with an accretion origin tend to be of lower metallicity and on highly elongated 
orbits with high energy. The GCs generally favor the extended outskirts of giant galax- 
ies (Lee et al. 1998, 2008b; Lee 2003; Dirsch et al. 2005), although they, on highly 
elongated orbits, must penetrate into the inner body of galaxies. This explains the 
observations that the bimodality in GC colors corresponds to a change in the kine- 
matic properties of the GCs, such that the blue GCs are dynamically hotter and there 
is a "bimodal" distribution of velocity dispersion. The demarcating radii between the 
inner, main part (the product of mergers) and the outskirts (the product of merger 
plus accretion) are not sharp and vary galaxy-to-galaxy depending on their individual 
histories of mergers and accretions. This is in line with the considerable diversity in 
the kinematics of the GC systems in giant elliptical galaxies (e.g., Hwang et al. 2008; 
Lee et al. 2010a). 
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Observationally, the intimately coupled histories of the GCs and constituent field stars 
of elliptical galaxies is encouraging in the context of mapping their metal contents, which is 
one of the most important yet least known properties of external galaxies. Even with the 
forthcoming lar ger telescopes, individual stars can be spatially resolved only in a few tens 
nearby galaxies ( lTolstoyll2006l ). Hence, we anticipate that the estimation of GC metallicities 
based on their colors, equipped with a proper correction for the GC-to-star offset, will become 
a preferred practical method in the quest to comprehend the evolution of galaxies beyond 
the nearby universe. Further refinement of the exact shape of the CMRs for GCs, and the 
resulting implications for the formation of GC systems and galaxies, should be the topic of 
much future work. 
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Fig. 1. — The color-magnitude diagrams, and eolor and inferred metallicity distributions for the 
GC system in giant elhptical galaxy, M87. (a) The g — z \s. u diagram for the GCs in M87. The 
HST/WFVC2 archival F336W (-u-band) images were used to derive u (ABMAG) for GC candi- 
dates. The -u-band catalog was matched with ACS/WFC g- and z-band (ABMAG) photometry of 
Jordan et al (2009) . The error bars represent the observational uncertainties of the corresponding 
magnitude bins. Jordan et al (2009) selected bona-fide GCs with their magnitudes, g — z colors, 
and sizes. We further employed color cuts in the u-band colors to filter out contaminating sources, 
especially background star-forming galaxies. The 591 GCs that have reliable u, g, and z measure- 
ments in common are only used, and the sample is n-band limited (uu < 0.2). (6) The HST/ ACS 
g — z color distribution for the M87 GCs. Distribution is normalized to the maximum value. Solid 
line is a smoothed histogram with a Gaussian kernel with a{g — z) = 0.05. (c) The relationship 
between g — z and [Fe/H] for the 40 low-extinction Galactic GCs (open circles), 33 M49 and M87 
GCs with ACSVCS photometry (blue circles), and 22 M49 and M87 GCs with SDSS photometry 
(red triangles). The references to the observed data used in the relation are summarized in Table 
1. The 5th-order polynomial fits to our model prediction for 13.9-Gyr GCs is overlaid (thick solid 
line). The a-element enhancement parameter, [a/Fe], is assumed to be 0.3. Thin grey straight line 
is for the linear least-squares fit to the data, {d, e, and /) The same as (ra c), but for the u — g 
colors. In (e), solid line is a smoothed histogram with a Gaussian kernel with cr(n — g) = 0.15. 
In (/), open circles, blue and red filled squares represent GCs in the Milky Way, M87, and NGC 
5128, respectively. The u — g colors of the GCs in the Milky Way and NGC 5128 were obtained 
from their U — B colors via the equation, {u — g) = 1.014 {U — B) + 1.372, derived from model 
data for synthetic GCs with combinations of age (10 ~ 15 Gyr of 1 Gyr intervals) and [Fe/H] (—2.5 
~ 0.5 dcx of 0.1 dex intervals), {g) The inferred GC MDFs based on the simple linear fit to the 
CMRs shown in (c) and (/). The red and blue hashed histograms are obtained using the g — z 
distribution in (6) and the v, — g distribution in (e), respectively. DistribTitions are normalized to 
the maximum values, (h) The same as (g), but with the inferred GC MDFs using the nonlinear 
color-to-metallicity transformations predicted by the YEPS models in (c) for the g — z distribution 
and in (/) for the u — g distribution. 
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Fig. 2.— The same as Figure 1, but for the M84 GCs. Their u (ABMAG) were obtained 
from HST/WFC3 archival images. The 306 GCs that have rehable u, g, and z measurements 
in common are only used, and the sample is ti-band limited (o"u < 0.2). In (6) and (e), the 
5th-order polynomial fits to our model prediction for 13-Gyr GCs is overlaid (thick solid 
line). Note that, in (/) the modeled u — g CMR differs from that in Figure If due to the 
difference in filter throughput between F336W's on WFPC2 and WFC. The main difference 
lies in that the WFG F336W has much less red leak than the WFPG2 F336W. The w - ^ 
colors of the GGs in the Milky Way and NGG 5128 were converted from their U — B colors 
via the equation, (u—g) — 1.296 (U — B) + 1.412, derived from model data for synthetic GGs 
with combinations of age (10 ~ 15 Gyr of 1 Gyr intervals) and [Fe/H] (—2.5 ~ 0.5 dex of 
0.1 dex intervals). The references to the observed data used in the relation are summarized 
in Table 1. 
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Fig. 3. — Color and metallicity distributions of GCs in ACSVCS galaxies, (a) The observed 
color histograms of GC sy stems for seven bi ns of host galaxy magnitude. The data are 
the same as in Figure 6 of iPeng et al.l (|2006[ ) and now listed in Table 4. The histograms 
are normalized by the GC number at their blue peaks, and multiplied by constants, C, for 
clarity. The magnitude bins are 1 mag wide and extend from Mb — —21.5 (—22 < Mb < 
-21, red, C = 1.0) to ~ -15.5 (-16 < Mb < -15, purple, C = 0.4). A Gaussian kernel 
of a{g — z) = 0.05 is applied, {h) Similar to Figures Ic and 2c, but the model predictions 
for various ages ranging from 10.5 Gyr (purple) to 13.5 (red) by equal age intervals of 0.5 
Gyr (see Table 2). The a-element enhancement parameter, [a/Fe], is assumed to be 0.3. 
The straight dotted line, again, is for the linear fit to the observational data points shown 
in Figures Ic and 2c, which are not plotted here for clarity, (c) The inferred MDFs of GCs 
in the same mag bins as in (a), multiplied by corresponding constants, C . Each GC MDF is 
obtained from the corresponding color distribution using the g — z io [Fe/H] conversion, and 
the same color code as in (a) is applied. The modeled ages for the host mag bins are given in 
the insets. In (c), the age is assume to be constant at 13 Gyr regardless of host luminosities 
between Mb — —15.5 (purple) and —21.5 (red). The grey dotted histogram represents the 
inferred MDF for the brightest (i.e., Mb — —21.5) bin based on the simple straight fit in 
(6), and is supposed to be compared to the red solid histogram, (d, e, /) The same as (c), 
but with age differences of 1.0, 2.0, and 3.0 Gyr, respectively, between the faintest {Mb — 
— 15.5, purple) and brightest (—21.5, red) bins. The brightest {Mb — —21.5, red) bin is set 
to be 13.5 Gyr. The data for the histograms in {d) are given in Table 5. 
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Fig. 4. — The same as Figure 3, but under a different assumption on tfie fiost age sequence 
from the brightest bin to faintest. It is assumed that the center bin, i.e., the fourth brightest 
{Mb — —18.5, green) bin, is set to be 13 Gyr. In (b), the YEPS model OMRs are displayed 
with ages ranging from 11.5 Gyr (purple) to 14.5 (red) by equal age intervals of 0.5 Gyr. 
The data for the histograms in (d), which is identical to Figure 3d, are listed in Table 5. 
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Fig. 5. — Comparison between MDFs of GCs and field stars. The inferred MDFs of ACSVCS 
GCs (thick curves) are taken from Figure 3d (see also Figure 4c? and Table 5), whereas 
those of field stars (grey histograms) are obtained from photometric studies on resolved 
red-giant stars in individual nearby galaxies. For a comparison, the dotted curve in each 
panel represents the MDF of ACSVCS GCs for the corresponding host luminosity bin that 
is derived based on the simple, straight fit shown in Figures Ic and 2c. The field star MDFs 
are measured in terms of [m/H] = Log(Z/Z0) - Log(X/X0), from which [Fe/H] is obtained 
using the equation, [Fe/ H] = LogfZ/Zr;,) - Log(X/X0) - 0.723 [a/Fe] = [m/H] - . 217, f or 
[a/Fe] = 0.3 (Table 2 of iKim et al.l jjoj, and s ee also Ishetrone. Cote, fc SargentI fcoOlh l 
The stellar metallicity is in the Zinn-West scale ( iZinn fc WestI 119841 ). All distributions a re 
normalized to the maximum values, (a) The stellar MDF of NGC 147 ( iHarris et al.ll2007a] Jbl) 
(a Local Group dwarf elliptical, Mb = —14.8) is compared to that of ACSVCS GCs in host 
galaxies with Mb — —15.5 (purple). (6) The sa me as (a), but for ho st galaxies in Mb — —19.5 
bin (yellow) and for field stars in NGC 3377 ( iHarris et al.ll2007al Jbl) (a Leo group elliptical, 
Mb = —19.2). The hashed area at [Fe/H] > —0.4 marks the region of low credibility on the 
stellar metallicity measurement due to photometric incompleteness, (c) The sar ne as (b), but 
for ho st galaxies in Mb — —20.5 bin (orange) and for field stars in NGC 3379 (IHarris et al. 
2007313) ( a Leo group elliptical, Mb = —20.6). {d) The same as (a), b ut for host galaxies i n 
Mb ~ -20.5 bin (orange) and for field stars in the halo of NGC 5128 jReikuba et al.ll2005[ ). 
a nearby field giant elliptical galaxy, with Mb = —20.9. The observed regions are ~ 21 and 
~ 31 Kpc away from the galaxy center, (e) The same as (d), but fo r host galaxies in Mr 
~ —21.5 bin (red) and for field stars in the halo of, again, NGC 5128 ( iRejkuba et al.ll2005l ). 
The observed region is ~ 8 Kpc away from the galaxy center. (/) The sar ne as (e), but for 
host galaxies in Mb — —21.5 bin (red) and for field stars in the halo of M87 ( iBird et al.ll2010l ) 
with Mb = —21.4. The observed regions are ~ 10 Kpc away from the galaxy center. The 
hashed area at [Fe/H] > —0.5 marks the region of low credibility on the stellar metallicity 
measurement due to photometric incompleteness. 
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Fig. 6. — Comparison between the MDFs of field stars of galaxies and their own GC systems 
for M87 (the top row), NGC 5128 (the second row), NGC 3377 (the third row), and NGC 
3379 (the bottom row). Among the galaxies in Figure 5, the four are those for which both 
the GC MDF and stellar MDF are currently available. The stellar MDFs (grey histograms) 
are identical to those in Figures 56, 5c, 5d, and 5/. The left-hand panels present the inferred 
GC MDFs based on the traditional linear color-to-metallicity transformations. By contrast, 
the right-hand panels show the inferred GC MDFs based on the inflected relations from the 
YEPS model. All distrib utions are riorma lized to the maximum values, (a) The stellar MDF 
(grey histogram) of M87 (IBird et al.ll2010l ) is compared to the MDF of 1745 GCs in the same 
galaxy (open thick histogram) obtained using the simple linear [Fe/H] vs. g — z relation in 
Figure Ic. The stars were sampled in an inner region (R ~ 10 Kpc) and the GCs were in the 
i nner halo (< 10 Kpc). (b) The stellar MDF (grey histogram) of NGC 5128 (jRejkuba et al. 



20051 ) is compared to the GC MDF (open thick hist ogram) obtained from C 



B — I histogram for 210 GCs in the same galaxy (IPeng. Ford, fc Freeman 



10 Bl anco 4-m 
2004a| M), usmg 



an empirical linear relation, [Fe/H] = 1.88 {B — I)q -4.11 (jSpitler. Forbes, fc Beasely 2008). 
The stars were sampled in an inner region (R ^ 20 and 30 Kpc) and the GC s were in the 
entire halo of this galaxy, (c) The stellar MDF (grey histogram) of NGC 3377 (IHarris et al. 



2007al lb[) is compared to the GC MDF (open thi ck histogram) o btained from HST/ ACS 



g — z histogram for 157 GCs in the same galaxy (jCho et al.ll201ll ) using the simple linear 
[Fe/H] vs. q — z rel ation in Figure Ic. {d) The stellar MDF (grey histogram) of NGC 3379 
faarris et al.ll2007al lbl) is compared to the GC MDF (open thick histogram) obtained from 
Rhode & Zepf (2004) 's B — R h istogram for 36 GC s with an improved empirical relation, 
[Fe/H] = 3.13 {B - R)o - 5.04 jHarris et al.ll2007bh . (e, /, g, and h) The same as {a-d), 
but for the GC MDFs of M87, NGC 5128, NGC 3377, and NGC 3379 are inferred from the 
YEPS model CMRs for g - z (13.5 Gyr), B - I (12.5 Gyr), g - z (12.8 Gyr), and B - R 
(12.5 Gyr), respectively. 
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Fig. 7. — Comparison between inferred GC MDFs and those from chemical enrichment 
models. The MDFs of ACSVCS GCs are taken from Figure 3d (see also Figure Ad and Table 
5). All distributions are normalized to the maximum values. (Upper) The MDF of ACSVCS 
GCs in the second brightest host bin with Mb — —20.5 (orange) is compared to that of a 
simple closed-box model with yield [Fe/H] = —0.69 (black solid line). For a comparison, 
dotted histogram represents the inferred MDF for the same Mb — —20.5 bin based on the 
simple straight fit shown in Figures Ic and 2c. (Lower) The same as the upper panel, but 
for the brightest host bin with Mb — —21.5 (red) and for model with yield [Fe/H] = —0.51 
(black solid line). For a comparison, dotted histogram represents the inferred MDF for the 
same Mb — —21.5 bin based on the simple straight fit shown in Figures Ic and 2c. 
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Table 1: References to the observational data for the [Fe/H] vs. color relations shown in 
Figures 1 and 2. 



Relations 


Galaxy Name 


References and Selection Criteria 








Spectroscopic [Fe/H] 


Broadband Color 




Figures Ic & 2c 


Milky Way 


1, 2 


1, 2 




The [Fe/H] vs. g — z relation 


M49 


1, 2 


1, 2 






M87 


1, 2 


1, 2 




Figure 1/ 


Milky Way 


3 


3 [{U -BY , E(B- 


V) < 0.3] 


The [Fe/H] vs. u — g relation 


NGC 5128 


4, 5 [i > 8 Gyr, S/N > 10] 


6 [{U -Bf] 




(w = HST/WFPC2 F336VF) 


M87 


1, 2 


This study [{u — g)] 




Figure 2/ 


Milky Way 


3 


3 [{U -Bf , E(B- 


V) < 0.3] 


The [Fe/H] vs. u — g relation 


NGC 5128 


4, 5 [t > 8 Gyr, S/N > 10] 


6 [(tZ-B)''] 




{u = HST/WFC3 F336VK) 


M87 


1, 2 


This study [{u — g)] 





"The equation, {u - g) = 1.014 {U - B) + 1.372, is used for _ff5'r/WFPC2 u-band F336W (see Figure 1). 
''The equation, (u - g) ^ 1.296 {U - B) + 1.412, is used for HST/WFCS w-band F336W (see Figure 2). 



References — ( 1) iPeng et al.l (120061): (2) Paper I; (3) Harris et a l (199 6 the 2010 edition); (4) 
Beaslev et all (|2008l ): (Sl lChung et al.l (|201ll ): (6) IPeng. Ford, fc FreemanI |2004alltf ). 
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Table 2: The g — z color from the YEPS model for GCs with ages (t) ranging from 9 Gyr to 
14 Gyr by steps of 0.5 Gyr. 



[Fe/H] 








9- z 












t = 9.0 


9.5 


10.0 


10.5 


11.0 11.5 


12.0 


12.5 


13.0 


13.5 


14.0 



-2 


5 





614 





623 





631 





652 





668 





696 





723 





745 
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Note. — The spectro-photometric model data (various colors and absorption indices) of the entire param- 
eter space are available at http:/7web.yonsei.ac.kr/cosmic7data/YEPSlhtm 
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Table 3: The median photometric errors in g — z and u — g oi five magnitude bins for the 
M87 and M84 GCs 



Mag bins 


Number of GCs 


g — z error 


u — g error 




M87 (m = 


= WFPC2 F336W) 




uo < 23.0 


20 




0.019 


0.024 


23.0 < Mo < 24.0 


87 




0.022 


0.028 


24.0 < Uo < 25.0 


172 




0.030 


0.054 


25.0 < Uo < 26.0 


247 




0.043 


0.110 


Uo > 26.0 


65 




0.058 


0.175 


Entire Sample 


591 




0.035 


0.078 




M84 (m 


= WFC3 F336W) 




Mo < 23.5 


9 




0.019 


0.021 


23.5 < Mo < 24.5 


65 




0.025 


0.032 


24.5 < Mo < 25.5 


107 




0.035 


0.059 


25.5 < Uo < 26.5 


109 




0.057 


0.111 


Mo > 26.5 


16 




0.082 


0.186 


Entire Sample 


306 




0.039 


0.068 
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Table 4: The g — z color histograms (Figures 3a & 4o) of GCs in ACS VCS galaxies binned 
by host galaxy 5-band magnitudes (Mb). 
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18 ± 5 
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± 
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19 ± 6 


12 ± 
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12 ± 5 
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9 


126 ± 11 
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± 
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± 
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4 ± 4 
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± 
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± 4 


± 3 
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± 
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1 ± 
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± 
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± 3 


-4 ± 3 


1.85 
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-3 ± 1 


-5 
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-7 ± 2 


-4 ± 
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-6 it 2 


-6 it 2 
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Table 5: The inferred MDFs (Figures 3d & 4rf, and Figures 5 & 7) of GCs in ACS VCS 
galaxies binned by host galaxy B-band magnitude {Mb)- 
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0.000 
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0.000 
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